Signatures of shape transition in odd-A neutron-rich Rubidium isotopes 
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The isotopic evolution of the ground-state nuclear shapes and the systematics of one-quasiproton 
configurations are studied in odd-A Rubidium isotopes. We use a selfconsistent Hartree-Fock- 
Bogoliubov formalism based on the Gogny energy density functional with two parametrizations, 
DIS and DIM, and implemented with the equal filling approximation. We find clear signatures 
of a sharp shape transition at A*" = 60 in both charge radii and spin-parity of the ground states, 
which are robust, consistent to each other, and in agreement with experiment. We point out that 
the combined analysis of these two observables could be used to predict unambiguously new regions 
where shape transitions might develop. 
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The study of the nuclear shape evolution as the number 
of nucleons changes is nowadays a highly topical issue in 
nuclear physics from both theoretical and experimental 
points of view (see for example Refs. and refer- 

ences therein). Specially interesting are those situations 
where the nuclear structure suffers drastic changes be- 
tween neighbor nuclides. These structural variations lead 
often to sudden changes of particular nuclear properties 
that can be used as signatures of phase/shape transi- 
tions 0, 0| . This is the case of the neutron-rich isotopes 
with masses A ~ 100. Intense experimental and 
theoretical [^, [l^, [lH efforts are being done to bet- 
ter characterize the structural evolution of the ground 
and excited states in this region of the nuclear chart. In 
particular, we have recently studied in Refs. 0, Q the 
structural evolution in even and odd neutron-rich Sr, Zr, 
and Mo isotopes. In this work we concentrate on the 
study of neutron-rich odd-mass Rb isotopes that have 
received considerable attention recently |12h15| and that 
exemplify a general pattern where structural fluctuations 
lead to observable effects. 

Being an odd-Z nucleus, the spin and parity of odd-A 
Rb isotopes are determined by the state occupied by the 
unpaired proton. The spectroscopic properties of the var- 
ious isotopes will be determined by the one-quasiproton 
configurations that, in principle, are expected to be 
rather stable against variations in the number of even 
neutrons. However, as it is known in neighboring nuclei, 
approaching TV ~ 60, the isotopes become well deformed 
[l|, |5| and the abrupt change in deformation induces sig- 
natures in nuclear bulk properties like the two-neutron 
separation energies and the nuclear charge radii, as well 
as in spectroscopic properties. In particular, the spin 
and parity of the nuclear ground state might change sud- 
denly from one isotope to another, reflecting the struc- 
tural change. 

We analyze in this work the bulk and spectroscopic 
properties of neutron-rich Rb odd isotopes within the 
selfconsistent Hartree-Fock-Bogoliubov (HFB) approxi- 
mation based on the finite range and density dependent 
Gogny energy density functional (EDF) fl6|. Previous 



studies in this region including triaxiality [s'] suggest that 
the 7-degree of freedom would not play a significant role 
in these isotopes and therefore, axial symmetry is as- 
sumed as a selfconsistent symmetry in this study. 

In addition to the well known DIS parametrization 
p7| of the Gogny-EDF, we also consider the most re- 
cent parametrization DIM fl8|. From this comparison 
we evaluate not only the robustness of our results, but 
we also explore the capability of DIM to account for the 
phenomenology of odd-A nuclei, not so well studied yet. 
The description of the odd-A nuclei involves additional 
difficulties because the blocking procedure requires the 
breaking of time-reversal invariance, making the calcula- 
tions arduous [l^, H^l • In the present study we use the 
equal filling approximation (EFA) , a prescription widely 
used in mean-field calculations to preserve the advan- 
tages of time-reversal invariance. In this approximation 
the odd nucleon sits half into a given orbital and half into 
its time-reversed partner. The microscopic justification 
of the EFA is based on standard ideas of quantum statis- 
tical mechanics 21} ■ The predictions arising from various 
treatments of the blocking have been studied in Ref. , 
concluding that the EFA is sufficiently precise for most 
practical applications. More details of our procedure can 
be found in Ref. 6]. 

The proton and neutron single-particle energies (SPEs) 
are shown in Fig. [T]as a function of the axial quadrupole 
moment Q20 for the even-even ^^Sr {Z — 38, N = 60) 
isotope. Fermi levels are plotted with thick dashed 
(red) lines. Asymptotic (Nilsson) quantum numbers 
[N,nz, A]K'^ are also shown in the proton case for the 
Q20 values where the energy minima are located (verti- 
cal arrows) in both oblate and prolate sectors. As one 
can see, the valence protons occupy the = 3 shell and, 
if deformed, they start to fill the 59/2 orbitals coming 
down from the iV = 4 shell. Neutrons occupy states be- 
longing to the iV = 4 shell and approaching iV = 60 
they start to populate the /iii/2 intruder orbitals com- 
ing down from the = 5 oscillator shell at large de- 
formations {Q20 4 — 5 b). These neutrons polarize 
the protons and drive them to populate the I59/2 shell 
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FIG. 2: (Color online) Experimental (a) excitation ener- 
gies and spin-parity assignments p^-ITsl. [23 ] compared with 
Gogny-DlS HFB-EFA results (b) for one-quasiproton states 
in odd-yl Rb isotopes. Prolate configurations are shown by 
black lines, oblate ones by red lines, and spherical ones by 
blue lines. 



FIG. f : (Color online) Single particle energies for protons and 
neutrons in '"'^Sr as a function of the axial quadrupole moment 
Q20. The Fermi level is depicted as a thick dashed red line. 
The results have been obtained with the Gogny-DlS EDF. 



at similar deformations. Thus, the underlying nuclear 
structure in this mass region is very sensitive to the oc- 
cupancy of these single-particle orbitals and the result 
is a rapid change in the nuclear deformation and in the 
spectroscopic properties as a function of both neutron 
and proton numbers. 

In particular, Rb isotopes approaching = 60 "will 
jump abruptly from slightly deformed oblate or prolate 
shapes to strongly prolate shapes. As a consequence, the 
unpaired proton that finally determines the spin and par- 
ity of the "whole nucleus in Rb isotopes, "will jump from 
the [303] 5/2" at low deformation to the [431] 3/2+ or- 
bital at large deformation. This is fully consistent "with 
the Jahn- Teller picture, "which states that nuclei avoid 
regions "with high single-particle level densities around 
the Fermi level, thus favoring deformation "where the en- 
ergy gaps are higher. One sees large energy gaps at large 



prolate deformations for the [431] 3/2+ orbital. This is 
also consistent "with the Federman-Pittel mechanism (2^ , 
"which points to the T = neutron-proton interaction as 
the driven force to"ward quadrupole deformation. This 
force is particularly intense bet"ween spin-orbit partners 
(neutron I57/2 and proton 1(79/2 in our case), as "well as 
bet"ween orbitals "with the same radial quantum numbers 
and large orbital angular momenta differing by one unit 
np = Tin and ip = ± 1 (neutron l/iii/2 and proton 
l(?g/2 in our case). 

In Fig. [21'we can see the experimental excitation ener- 
gies and spin-parity assignments (a) in odd-A neutron- 
rich Rb isotopes. They are compared to the one- 
quasiproton states predicted by the Gogny-DlS HFB- 
EFA calculation (b). The excited states in a given isotope 
are referred to the corresponding ground state, regard- 
less its shape. Prolate configurations in our calculations 
are sho"wn by black lines, oblate ones by red lines, and 
spherical ones by blue lines. The quasiparticle states are 
labeled by their K'^ quantum numbers. The most impor- 
tant configurations arc joined by dashed lines follo"wing 
the isotopic evolution. In addition, the ground states are 
labeled by their asymptotic quantum numbers. 
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FIG. 3: (Color online) Gogny-DlS excitation energies of 
single-quasiproton prolate (a) and oblate (b) states in Rb iso- 
topes. Hole states are plotted below zero energy and particle 
states are plotted above. The absolute ground states are in- 
dicated with a circle. 
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FIG. 4: (Color online) Same as in Fig. El but for Gogny DIM. 



Experimentally [T3-[T1 [H one observes ,r ^ 3/2" 
ground states in ^'^'^^'^^Rb, then 5/2" states in ^^-SSj^^^ 
and finally 3/2+ states in the heavier isotopes. The 9/2+ 
states appear all the way up to ^^Rb as excited states. 
The most striking feature observed is a jump at iV = 
60 {A = 97) from 5/2" to 3/2+ ground states. 
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FIG. 5: (Color online) Isotopic evolution of the quadrupole 
deformation parameter /3 of the energy minima obtained from 
Gogny DIS and DIM calculations. 



The theoretical interpretation of these features can be 
understood from the analysis of our results in the lower 
panel (b). Rb isotopes evolve from spherical shapes in 
85-89 around TV = 50 with the spherical /5/2, P3/2, and 
Pi/2 shells involved, to slightly deformed shapes (oblate 
and prolate almost degenerate) in ^-'^"^^Rb, and finally to 
well deformed prolate shapes in ^^"^"^Rb. In the lighter 
isotopes the three spherical shells mentioned above are 
very close in energy and thus, we can see all the split K'^ 
levels coming from them at very close low-lying excita- 
tion energies. In the case of ^^~^^Rb isotopes we obtain 
oblate ground states with both oblate and prolate low- 
lying excited states very close in energy. For instance, we 
see that the experimental 3/2~ ground states in *^'^^Rb 
appear in the calculation as excited states at less than 
0.1 McV. In the case of ^^'^^Rb the experimental ground 
states 5/2~ appear as prolate excited states, while the 
calculations produce oblate ground states. The ground 
states of the heavier isotopes are predicted to be 3/2+ in 
agreement with experiment. The calculations also pre- 
dict excited low- lying 9/2+ states with an oblate nature, 
which are observed experimentally. 

In the next figures we compare the spectroscopic prop- 
erties of the Gogny DIS and DIM parametrizations. 
Contrary to Fig. in Fig. |3]for DIS and in Fig. |3]for 
DIM, we have separated the prolate (a) and oblate (b) 
states and have plotted the most relevant hole states be- 
low zero energy and the particle states above. The abso- 
lute ground states, either oblate or prolate, are indicated 
with a circle. Specifically, in the figures corresponding to 
the prolate shape, we have plotted the evolution of the 
3/2+((79/2) states, which are ground states for ^''~-^°^Rb. 
In the upper region we find the 5/2+(gg/2), while in 
the lower region we have the 1/2+ (59/2), 3/2^(^3/2), 
l/2^(/5/2), and 7/2~ coming up from the f'j/2 shell. In 
the oblate case (b) we can see the 9/2+ (179/2) state, which 
is ground state in ^^Rb, the 7/2+((jr9/2) as particle state 
and the l/2~ and 3/2~ from /5/2 as hole states. 

Very similar results are obtained from both 
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parametrizations, which answer our original ques- 
tions about robustness of the calculations and reliability 
of DIM. The only difference worth mentioning is that 
DIM produces slightly lower excited states. This feature 
can be understood from its larger effective mass that 
makes the single particle spectrum somewhat more 
dense with DIM. 

To further illustrate the role of deformation and spin- 
parity assignments in the isotopic evolution, we display 
in Fig. [5] the axial quadrupole deformation /3 of the en- 
ergy minima as a function of N with both DIS and DIM 
parametrizations. The deformation of the ground state 
for each isotope is encircled. We can see that starting 
at the semi-magic isotope with = 50 we get two min- 
ima in the prolate and oblate sectors. These are 5/2^ in 
the prolate case and l/2~ (or 9/2+ very close in energy) 
in the oblate case. The calculations with the two forces 
predict oblate ground states, but with prolate solutions 
which are very close in energy and that could perfectly 
become ground states for slightly different type of cal- 
culations. Thus, in ^^Rb the l/2~ and 5/2 ~ are practi- 
cally degenerate (see Fig. [2]). In ^^Rb (^'^Rb) they are 
separated by less than 0.2 (0.5) MeV. The experimental 
information seems to favor prolate solutions since 5/2 ~ 
ground states are observed. The ground states of heav- 
ier isotopes starting at ^^Rb are 3/2+ in agreement with 
experiment and strongly prolate. The 9/2+ oblate con- 
figurations appear all the way toward heavier isotopes. 
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FIG. 6: Gogny-DlS HFB results for d(rl) (a) and S2n (b) 
in odd-A Rb isotopes compared to ex per imental data from 
Ref. [2^, [2^ for masses and from Ref. [27l ] for radii. Results 
for prolate, oblate, and spherical minima are displayed with 
different symbols (see legend). Open circles correspond to 
ground-state results. 

Finally, we show in Fig. [5] for DIS and in Fig. 
[7] for DIM, the charge radii differences (a) defined as 
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FIG. 7: Same as in Fig. El but for Gogny DIM. 



rections as in Ref. [5], and the two-neutron separation 
energies S2n (b). We compare our results for S{rl) with 
isotope shifts from laser spectroscopy experiments 12711 
and our results for S2n with masses from Refs. [1^ |26| . 
The results for S2n agree in general with the measure- 
ments and are very similar between DIS and DIM. A 
somewhat better agreement is found with DIM at the 
beginning of the shell {N > 50). In any case, it would 
be interesting to see whether the incipient experimental 
tendency to raise the S2n values in the heavier isotopes 
and not predicted by the calculations, still persists in 
more exotic nuclei. In this respect, we should remark 
that for the heavier Rb isotopes the mass data are not 
completely reliable because they are based on /3-endpoint 
measurements. It has been shown that this type of indi- 
rect measurements systematically tend to underestimate 
the Q-values leading to very strong binding. This prob- 
lem has been discussed in Ref. [28| for neutron-rich Y 
and Nb isotopes. Thus, it would be very helpful to ex- 
tend and improve mass measurements reducing the still 
large uncertainties in neutron-rich Rb isotopes. 

The nuclear charge radii plotted in the upper figures 
show a jump at TV = 60, where the radius suddenly 
increases. This experimental observation is well repro- 
duced in our calculations, where the encircled ground 
states show that the jump occurs between iV = 58 and 
iV = 60. This jump correlates well with the one observed 
in the spin-parity of the ground states displayed in Fig. 
dfrom 5/2" to 3/2+. Theoretically, the jump in S{rl) 
could be understood as the result of the change from an 
oblate to a prolate shape. But it could also be inter- 
preted as a jump from a slightly prolate shape to a well 
deformed prolate nucleus, because the oblate and prolate 
shapes in the transitional nuclei A'^ = 54 — 58 have very 
similar absolute values of deformation (see Fig. [5]). This 
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similarity gives raise to practically the same charge radii, 
as it can be seen by comparing the prolate (upward red 
triangles) and oblate (downward blue triangles) radii in 
this region. The only difference worth mentioning be- 
tween DIS and DIM is the oblate radii at = 60 which 
corresponds to different configurations (7/2+ in DIS and 
9/2+ in DIM) and different deformations (see Fig. [5])- 

Then, we see that contrary to S2m which does not sig- 
nal any clear signature in Rb isotopes, S{r'^) clearly in- 
dicates the existence of a shape transition in agreement 
with experiment and well correlated to the spin-parity 
jump. It would be also very interesting to extend the iso- 
tope shifts measurements, which are rather old, to more 
exotic isotopes to confirm also the prolate stabilization. 

In summary, we have studied the shape evolution 
in odd-A Rb isotopes from microscopic selfconsistent 
Gogny-EDF HFB-EFA calculations. We have analyzed 
various sensitive nuclear observables, such as two-neutron 
separation energies, charge radii, and the spin-parity of 
ground states in a search for signatures of shape tran- 
sitions. We have found that, although the masses are 
not very sensitive to these shape changes in Rb isotopes, 
the charge radii and the spin-parity are. In addition, 
the signatures found are all consistent to each other and 
point unambiguously to the existence of a shape transi- 
tion. The correlations between all of these signatures are 
useful to establish a protocol to look for shape transitions 
in the future in other regions of the nuclear chart and to 
predict and identify them. We have also analyzed the iso- 
topic evolution of the one-quasiproton configurations and 
have compared the predictions of two different Gogny- 
EDF parametrizations, demonstrating the robustness of 
our calculations. 

The experimental information available in neutron-rich 
Rb isotopes and in general in this mass region is still 
very limited. It would be highly desirable to extend the 
experimental programs for mass, charge radii, and spec- 
troscopic measurements to these exotic regions at exist- 
ing facilities like ISOLTRAP 29] at ISOLDE/CERN and 
IGISOL jsoj at the University of Jyvaskyla, or at fu- 
ture ones like MATS and LaSpec Q at the FAIR facility, 
where we can learn much about structural evolution in 
nuclear systems. 

The theoretical efforts should be also pushed forward 
by improving the formalism including triaxial degrees of 
freedom in those regions where this can be an issue or 
dealing with the odd systems with exact blocking treat- 
ments. The quality of our mean field description could 
be improved from configuration mixing calculations in 
the spirit of the Generator Coordinate Method with the 
quadrupole moment as generator coordinate. Neverthe- 
less, such an approach for odd nuclei is out of the scope 
of the present study. It is already known ^] that such a 
configuration mixing reduces the jump in S2n predicted 
in pure mean field approximations when crossing shell 
closures, improving the agreement with experiment. This 
could be particularly relevant for the light isotopes con- 
sidered in the present study, where the spherical minima 



are rather shallow. For heavier isotopes, the two min- 
ima, oblate and prolate, are separated by spherical bar- 
riers of about 3 MeV and appear about 1 MeV apart. 
The effect here is not expected to be significant because 
a single shape would be enough to account for the prop- 
erties studied. A simple two-state mixing model could 
be used to mix the intrinsic deformed configurations into 
the physical states [Ij. Assuming that the wave functions 
are localized at largely different values of the deformation 
coordinate, one can neglect the cross terms. The final ef- 
fect would be to place the value of the radius in between 
the oblate and prolate values according to the mixing 
amplitude of those intrinsic configurations into the phys- 
ical wave function. The consequence would be a slightly 
smoother transition. 
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